Crossover and scaling in a nearly antiferromagnetic Fermi liquid in two dimensions 
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We consider two-dimensional Fermi liquids in the vicinity of a quantum transition to a phase with 
commensurate, antiferromagnetic long-range order. Depending upon the Fermi surface topology, 
mean-field spin-density-wave theory predicts two different types of such transitions, with mean- 
field dynamic critical exponents z = 1 (when the Fermi surface does not cross the magnetic zone 
boundary, type A) and z = 2 (when the Fermi surface crosses the magnetic zone boundary, type 
B). The type A system only displays z = 1 behavior at all energies and its scaling properties 
are similar (though not identical) to those of an insulating Heisenberg antiferromagnet. Under 
suitable conditions precisely stated in this paper, the type B system displays a crossover from 
relaxational behavior at low energies to type A behavior at high energies. A scaling hypothesis is 
proposed to describe this crossover: we postulate a universal scaling function which determines the 
entire, temperature-, wavevector-, and frequency-dependent, dynamic, staggered spin susceptibility 
in terms of 4 measurable, T — 0, parameters (determining the distance, energy, and order parameter 
scales, plus one crossover parameter). The scaling function contains the full scaling behavior in all 
regimes for both type A and B systems. The crossover behavior of the uniform susceptibility and 
the specific heat is somewhat more complicated and is also discussed. Explicit computation of the 
crossover functions is carried out in a large N expansion on a mean-field model. Some new results 
for the critical properties on the ordered side of the transition are also obtained in a spin-density 
wave formalism. The possible relevance of our results to the doped cuprate compounds is briefly 
discussed. 
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I. INTRODUCTION 

A number of recent works [|J^8| have proposed a description of the low temperature spin dynamics of strongly 
correlated electronic systems using their proximity to an actual or hypothetical zero temperature, magnetic, quantum 
phase transition. Such an approach has the advantage of allowing development of a systematic expansion about a 
point with non-trivial spin correlations and with strongly interacting excitations. Some of these proposals have arisen 
in the context of two-dimensional interacting electron models of the doped cuprate superconductors, while others 
considered the heavy- fermion compounds. The ideas of this paper will be presented in the former context, although 
our approach may be of a more general utility. 

We begin our discussion by reviewing two recent complementary approaches to magnetic quantum transitions in 
the cuprate compounds. A unified picture of spin fluctuations in the lightly-doped cuprates obtained from these past 
approaches, and our present work, is provided towards the end of this paper in Sec. M. We wish to emphasize at the 
outset that all of the discussion in this paper refers to spin fluctuations associated with commensurate antiferromag- 
netic ordering; for the cuprates this corresponds to an ordering wavevector Q — (w,iv). We have little to say here 
about the case of incommensurate ordering. 

The first approach begins from the insulating parent compound, like LCZ2CMO4, whose spin fluctuations can be 
modeled by a spin-1/2 Heisenberg antiferromagnet on a square lattice JjJ. The low energy excitations of the antifer- 
romagnet are believed to be well-described by a continuum 0(3) non- linear sigma model field theory. The non-linear 
sigma model is parameterized by a single coupling constant, g, which measures the strength of quantum fluctuations 
in the system; for g < g c , the system has Neel order, while for g > g c , it is in the quantum disordered phase. The 
transition at g = g c has been studied in some detail [|4J — it has dynamic critical exponent z = 1 and leads to a 
quantum disordered phase in which at T = the low energy magnon excitations have a gap and an infinite lifetime. 
The current experimental and theoretical consensus is that the S = 1/2 Heisenberg antiferromagnet has macroscopic 
Neel order in the ground state, and therefore should map onto a sigma model with g < g c . Now consider doping 
this antiferromagnet with a small number of holes. At very small doping, the holes form small elliptical pockets 
at (ir/2,ir/2) and symmetry related points in the magnetic Brillouin zone p-[f l| . Quantum fluctuations associated 
with these holes decrease the spin stiffness of the antiferromagnet J9) and should therefore drive the effective value 
of g closer to the quantum transition point at g = g c . (We are assuming here that conditions are such that there 
is a direct transition from a commensurate long-range ordered state to a commensurate quantum-disordered state 
with increasing doping E3] ; we are neglecting the possibility of an intermediate incommensurate long-range ordered 
state |13|,|l4| ) • Suppose that quantum fluctuations are strong enough that the disordering transition occurs while the 
holes still occupy small pockets. The critical properties of such a transition were studied by one of us 113] in the 



framework of the Shraiman-Siggia model 13 1: it was found that the dynamic critical exponent remained at z = 1, and 
the mobile charge carriers only introduced a small, but negligible damping of the magnon excitations at the critical 
point. The full structure of the quantum disordered phase in this model, and in particular, the topology of its Fermi 



surface is not well understood: we will discuss these issues further in Sec [II C. 

The second approach begins with the opposite limit of large doping, where the dilute system of electrons is pre- 
sumably well described as a Fermi liquid. Using mean-field spin-density-wave ideas, a scenario for the onset of 
antiferromagnetic order in such a Fermi liquid was proposed many years ago by Hertz fll5fl , and extended recently by 
Millis H . It turns out to be important to distinguish two cases [pf depending upon the value the ordering wavevector 
Q (for our case Q = (it, 7t)), and the shape of the Fermi surface in the quantum-disordered phase. 
(A) Damping of spin excitations with momenta near Q, due to conversion into particle-hole pairs, is forbidden: This 
is the case when Q cannot connect two points on the Fermi surface. For a circular Fermi surface this corresponds to 
Q > 2kp where kp is the Fermi wavevector. The transition then has the mean-field exponent 2 = 1 and is rather 
similar to the discussion above on the first approach. The scaling results of Ref. Q apply most ly unchanged: some 



simple modifications are necessary for the uniform susceptibility, and are discussed in Sec. 1VA 
(B) Damping of spin excitations with momenta near Q, due to conversion into particle-hole pairs, is allowed: This 
is the case when Q can connect two points on the Fermi surface, which for a circular Fermi surface corresponds to 
Q < 2kp. This transition has z = 2 in mean-field theory. Accordingly, the magnon excitations in the quantum 
disordered phase are overdamped and relaxational. 

We will not discuss the special case Q = 2kp in this paper. 

A phenomenological form for the magnetic susceptibility near a type B transition was introduced in context of 
cuprate superconductors by Millis, Monien and Pines |16| to explain NMR data in YBa2Cu^O-j. The z = 1 to z = 2 
crossover with doping was later proposed by Pines and one of us H and by Barzykin et al. M to describe evolution 
of the experimental NMR and neutron scattering data with oxygen concentration for YBo2Cu^Oq+ x . Very recently, 



Liu and Su || considered a two-component "Kondo-lattice"-type model of the cuprates, with separate localized spin 
and itinerant electron degrees of freedom. They assumed conditions that were appropriate to have a type B phase 
transition with z = 2, and found behavior characteristic of this value of z at the lowest energies or temperatures. At 
higher temperatures however, they found a crossover to behavior characteristic of z = 1. Earlier, the temperature- 
induced crossover to the z = 1 regime was discussed in Ref (I?]] in relation to YBa2Cu3Oe.es- It is not difficult to see 
that this crossover is in fact a rather general phenomenon for type B systems — there will always be reactive terms 
~ lu 2 (to is a measuring frequency) in spin response functions which will overwhelm dissipative terms at large enough 
lu. The crossover should be especially pronounced in systems where the damping constant, 7, is small. 

The main purpose of this paper is to discuss the scaling properties of the quantum disordered phase of systems near 
transitions of type B. We will also make a few comments about type A systems to which the analysis in Ref. [0 will 
mostly apply. We will assume that conditions in the type B system are such that the higher energy crossover to type 
A behavior occurs at energies which are significantly smaller than other high energy cutoffs like the Fermi energy or 
the exchange constant. (For type B systems which violate this condition, our results reproduce the correct asymptotic 
critical singularities, but make a particular choice for numerical scale-factors which are in fact non-universal in this 
case; these statements will be made more precise later). We also assume, of course, that no other unrelated low energy 
scale appears. We will propose a scaling hypothesis, in which the wavevector (q), frequency (lu), and temperature (T) 
dependence of the staggered spin susceptibility can be expressed in terms of a scaling function which involves only 
four experimentally determinable, T — input parameters. Only these four parameters are dependent on the details 
of the microscopic interactions in the ground state; everything else is universal and can, in principle, be computed in 
a long distance field theory. The general scaling arguments will be presented in the Sec |l[ The crossover behavior 
of the uniform susceptibility and the specific heat is somewhat more complicated and requires additional microscopic 



parameters — this is discussed in Sees IV A and IV B 



We will then illustrate these scaling ideas in two model calculations. 



In Section [II we consider magnetic phase transitions in a standard spin-density-wave (SDW) formalism. We will 
examine a simple model — a Hubbard model with first (t) and second neighbor hopping (£') for the fermions, which 
displays a transition of either type A or type B, depending on the ratio t/t' . Note that in both cases, the Fermi surface 
in the quantum-disordered phase is large (i.e. encloses a volume given by the total number of electrons). The SDW 
results for the magnetic susceptibility will illustrate the important differences in the nature of the spin excitations 
between type A and B transitions. The SDW theory will also be used to obtain new mean-field results for critical 
properties on the magnetically ordered side of the transition. Finally, we will briefly discuss the relation between the 
SDW results and those for the fluctuation-driven magnetic transition within the Shraiman-Siggia model. 



In Section IV, we first use the results of Section III to motivate a model field theory, S, to describe the quantum- 
disordered phase in case B and its crossover to type A behavior. The model S, which relies on a mean-field ap- 
proximation for the fermions, will turn out to be precisely the one proposed by Liu and Su H; we caution however 
that the validity of S as a description of the underlying fermionic excitations has not been conclusively established 
asymptotically close to the critical point. We will use a large N expansion to compute explicit results for the crossover 
scaling functions of S. 

Finally in Section [V] we will discuss the experimental relevance of our results and state our main conclusions. 

II. SCALING HYPOTHESES 

We now present some general scaling ideas to describe the crossover between type A and type B transitions. In 
the terminology of the well-developed theory of crossover phenomena between two critical points in classical phase 
transitions 118UT9], we need to distinguish between the primary and secondary fixed points: the primary fixed point 
has 2 relevant directions, while the secondary point has only a single, relevant direction associated with the "thermal" 
operator which drives one across the transition. In our case, it is clear that the primary fixed point is the type A fixed 
point with z = 1. The mean-field value of z — 1 for type A is expected to be robust and not suffer any corrections 
from fluctuations, as the coupling between the spin-fluctuations and fermionic quasiparticles is quite weak in this case. 
Adding any damping to the magnon excitations should be a relevant perturbation at this fixed point as it introduces 
additional low energy excitations (it could be dangerously irrelevant, a possibility which we shall briefly refer to later, 
but not explore in any detail). The coupling, analogous to g, which tunes the system across the transition is the other 
relevant perturbation. The type B fixed point must therefore be secondary. (We are also assuming here that the effect 
of the damping is associated with only a single relevant operator; there could be more than one. It should be easy to 
extend the following scaling analysis to this case, but we shall not do it in the interests of simplicity). An important 



feature of the standard crossover theory 118 19 is that the crossover scaling functions between the two fixed points 



are expressed in terms of eigenoperators and exponents of the primary fixed point, while the critical singularities of 
the secondary fixed point appear as nonanalytic behavior in the crossover functions themselves. This result forms the 
basis of our analysis below. 

We will restrict our scaling results below to the to the quantum disordered side of the transition. We begin our 
discussion by introducing the four parameters which will characterize the quantum-disordered ground state; the finite 
temperature properties of the staggered spin susceptibility, \s , will then be described by universal scaling functions 
of these four parameters only. Imagine we have available, either through experiments or computer simulations, the 
T = value of the imaginary part of the local, on-site, dynamic spin susceptibility Xl( w ) °f the system of interest; 
the local susceptibility is obtained by integrating the dynamic susceptibility over momenta in the vicinity of Q (where 
it equals % s (g, w)). In an insulating antiferromagnet, whose transition is described completely by the primary, z = 1 
fixed point, x'L wou ld have the form shown in Fig nk B: there is gap below which the spectral density is strictly zero, 
a discontinuity at the gap, and for large u> we have x'L ~ ^ where r\ is an anomalous dimension of the z — \ fixed 
point [|3|j4|. Here, by large u, we mean frequencies which are large compared to the spin fluctuation energy scale, but 
small compared to upper cutoffs like the an exchange constant or the Fermi energy; this and similar restrictions will be 
implicitly assumed in the remainder of the paper. Let us now examine the change in the spectrum due to a relevant 
perturbation which moves the system towards the secondary, type B, fixed point. Under appropriate conditions, 
mobile fermionic carriers can act as such a perturbation, and the damping due to the particle-hole continuum will 
introduce some sub-gap absorption: x'L wm then look like Fig |l|b. The gap has turned into a pseudo-gap, and x'L ~ w 
for small ui. However, we will still have x'L ~ W?J f° r kirge w as the primary fixed point behavior is expected to 
continue to hold at large u>. We now extract three parameters from this form for x'L'- 
(i) An overall amplitude Z: This is defined by 

Z = 4 lim $M at T = 0. (2.1) 

The factor of 4 is for convenience in a later model calculation. 

(ii) and (Hi) Two energy scales, A and T, which measure the pseudogap and the strength of the damping respectively. 

These are obtained by solving the constraints 

x 'L(hu, = A) = ^ ; i- = J-?- lim 4^ bothatT = 0. (2.2) 



The factor of 8 in the first equation is chosen to be exactly twice the factor 4 in (2.1); the factor of 47r in the second 
equation is arbitrary and is for future convenience. The parameter M is the energy below which the anomalous 
scaling of the primary fixed point x'L ~ ^ stops, and a convenient choice is M = amax(r, A), where a is of order, 
but larger than, unity, the factors of M v above also ensure that A and F have units of energy. Because of the tiny 
value of r\ « 0.03, the M v can be dropped while determining T and A from data; this is however not true in any 
analytical computation of scaling functions, as they are essential in canceling cut-off dependencies and obtaining a 
universal result. The pseudo-gap A has therefore been defined as the frequency at which x'L ^ a ^ s to half its large 
frequency value (modulo factors of M v ), while T is determined from the sub-gap absorption. For T < A, A is roughly 
the location of the knee in x^(w), while T is its width (see Fig [j]b); however, the definitions above hold for all T/A. 
For very small T, when the system is very close to the primary fixed point, T and A measure the strength of the two 
relevant perturbations away from it; this is no longer true for large T/A, but our scaling results below will continue 
to be valid. 

(iv) The final parameter sets the normalization of length scales. A convenient choice is to use the T = spin wave 
velocity, c, defined by the q dependence of the peak in the imaginary part of the staggered spin susceptibility Xs (q, ui) 
at large q (In all expressions for Xs only, it is implicitly assumed that wavevectors, q, are deviations from the ordering 
wavevector Q). 

Scaling functions of the primary z = 1 fixed point were studied in some detail in Ref . B , and we can easily extend 
that analysis to write down the crossover scaling functions for the present case. The main difference is that we have 
a new relevant parameter, T, which must be included in the scaling functions. We consider only case of the scaling 
function of the staggered susceptibility, Xs, here, leaving for later other observables which have a somewhat more 
complicated crossover behavior. Using the fact that T was defined above to have the dimensions of energy, we can 
write 

he 



^^MPte) M^A,F) (2.3) 



where $ s is a fully universal, dimensionless crossover function, and the dimensionless arguments measure values of 
the parameters in units of fc^T: thus 

__ hcq _ _ huj y- _ A — _ r . . 

ksT ksT ksT ksT 

All letter exponents rj, v (to be used later) refer to the primary fixed point. The main condition for the validity of 
( |2.3D is that the five energy scales, Tick, Tilu, A, T and fc^T are all significantly smaller than upper cutoffs like the 
Fermi energy or an exchange constant. The ratios of the five energy scales can however be arbitrary. 

At r = 0, <& s reduces to the result of Ref H (there are some minor differences in conventions), and we obtain 
the scaling function for case A. The criticality of the secondary fixed point appears in the T — > oo limit. Then $ s 
should collapse into a secondary scaling function $f [|18|,|l9|] with the T argument removed. The other arguments, 
and the overall scale, of $ s will be multiplied by powers of V, so that <&~ and its arguments have scaling dimensions 
appropriate to the secondary fixed point. To illustrate this point more clearly, let us assume that the structure of 
the type B fixed point is identical to the Gaussian z = 2 fixed point obtained in spin-density-wave mean-field theory; 
then we will have 

lim $ s (fc,uJ, A,T) = _ $f I . uJ, —-— I upto logarithmic corrections. (2.5) 

r^oo r n Vr ' V ' J 

The logarithmic corrections arise because the z — 2 critical point has marginal perturbations in d = 2 pa]. It is 
however important to note that these logarithmic corrections, and indeed the entire leading critical behavior of the 
secondary fixed point, are contained completely within the primary universal scaling function $ s . Everything about 
the logarithmic terms is universal and, in particular, there is no non-universal cutoff-dependent argument to the 
logarithms. This universality is a key point, and is a consequence of our earlier choice to define the scaling functions 
with respect to the primary fixed point. In systems with T so large that it is greater than a high energy cutoff like 
the Fermi energy, this strong universality will not hold; however, the large T limit of $ s can still be used to describe 
the critical behavior, with the understanding that the various numerical scale factors are not correct as they are now 
non-universal. 

III. MEAN-FIELD SDW CALCULATIONS FOR THE HUBBARD MODEL 

The remainder of the paper will make the rather abstract arguments of Section |] concrete by presenting explicit 
calculations in some simple models. We start with the mean-field SDW theory of the disordering transition in doped 
antiferromagnets. The simplest model which displays such transition is a one-band Hubbard model given by 

Here j and j' label the nearest and the next-nearest neighbors, respectively, and n = c)c is the particle density. 
Depending on the density of carriers and the ratio t'/t, the Fermi surface of free electrons can either be closed, in 
which case it is centered at (0, 0) and located entirely inside the magnetic Brillouin zone, or it can be an open Fermi 
surface which crosses the magnetic Brillouin zone boundary (Fig. Q). We will not specify the values of t and t' for which 
the Fermi surface is open or closed for a particular doping concentration, but rather consider the critical behavior of 
the dynamic spin susceptibility in both cases. 

Let us assume that the model has commensurate antiferromagnetic long-range order down to a transition point 
(this assumption may not hold for a particular choice of t and t', but it should always hold for related models with 
additional momentum dependence in U ]20| ; this momentum dependence however does not lead to new physics near 
the transition, and we will not consider this complication here). The (jr, it) ordering implies that, e.g., the z component 
of the spin-density operator, S(q) — (1/2) ^ fc ajj + „ a (?a.,i3Q>k,f)i has a nonzero expectation value at q = Q = (tt, 7f). In 
the SDW approach Q, the relation (X^fc a l+Q T a fe,t) = — (J2k a k+o i a M) = (Sz), is use d to decouple the quartic 



term in (3.1). After decoupling and diagonalization of the quadratic form, one obtains 



H S vw = ^ K<L C *° + E tAo d ^ (3-2) 

k 

where primes on the summation signs indicate that it is over the reduced magnetic Brillouin zone, and we introduced 



E C k=E^+4, E d k = -E-+e+ (3.3) 

where E^ = JNq + (e^) 2 , and N = U(S Z ), e^ = (e fc ±e/c + Q)/2 , e k = -2t(cos k x + cos k y ) - 4t' cos k x cosk y . We will 
refer to the quasiparticles described by c and d operators as conduction and valence fermions, respectively. Finally, 
the self-consistency condition on (5*2) is 



d _ 

(3.4) 



X^' n k~ n k 



u r E * ' 

Below we will assume that t' is negative. This is consistent with numerical calculations g3|,Q and fits to the 
measured |2[| shape of Fermi surface for YBaiCu-sOi. 

We now describe the evolution of the Fermi surface with doping. The shape of the Fermi surface is determined from 
E c ' d = p. At half-filling, the system is an insulator (\fi\ < No), all valence band states are occupied and all conduction 
band states are empty. At finite doping, the chemical potential moves into the valence band. As the maximum of 
E d is located at (ir/2,ir/2) and symmetry related points, the Fermi surface first opens in the form of hole pockets 
located around these points (Fig Br,) . The subsequent evolution of of the Fermi surface with doping in the SDW phase 
depends on the ratio t'/t. For sufficiently large \t'\, the bottom of the conduction band also becomes smaller than 
\x above some particular doping concentration, and electron pockets appear in the conduction band (Fig pb). The 
energy gap between the hole and electron pockets shrinks gradually as No decreases, and, at the transition point, the 
electron and hole pockets merge at (xq, it— xq) and symmetry related points, where cos 2 xo = |/i/4i'|, and the value of 
|/i| is determined from the solution of Eq. (3.4) with Nq — * (Fig pp). It can be verified that this transition is of type 



B, and Q can span two points on the Fermi surface in the paramagnetic phase. For small enough \t'\, the situation is 
different: the size of hole pockets continues to increase with doping, and at some concentration x\, the hole pockets 
extend to (n, 0) and related points while magnetic order is still present (Fig pp). At concentrations smaller than x\, 
electron pockets either do not appear at all, or they first appear, their size passes through a maximum, and then 
they disappear at some doping concentration x 2 < X\. At some concentration X3 > X2, the long-range order finally 
disappears, but the Fermi surface at the transition is located inside magnetic Brillouin zone boundary (Fig He). The 
disordering transition is therefore of type A |2q] . 

We now present results from a mean-field SDW calculation of the uniform and staggered susceptibilities; details 
can be found in the appendix. We shall mainly be interested in the small q and u> behavior of these responses when 
No is small. We shall find that an important feature of type B systems is that the limits q, lo — ► and No — > do 
not commute. We consider the two cases separately. 

A. Type A 

The transverse, staggered, susceptibility is found to behave as 

iV 2 1 

Xs_L(q,u) = — - — 2—2 : -2 (3.5) 

p s g z — lo j c — laujqjyiQ 

where a, c are constants (i.e. finite as No — ► 0) and p s ~ Nq ~ (g c — g) is a spin stiffness. It is clear that the damping 
is negligible as Nq — > 0, and that the limits q, lo — » and A?o — ► do commute. The susceptibility is also that expected 
for a z — 1 system. 

In the quantum disordered phase we then have %s ~ [q 2 — lo 2 /c 2 + A 2 /c 2 ) -1 , where the gap A vanishes at the 
transition point. The pole in Xs corresponds to a positive energy, spin-1 particle- hole bound state which has center 
of mass momentum near Q, and which lies below the bottom of the particle- hole continuum at these momenta. 

The transverse, uniform susceptibility at No — > can be represented as a sum of the interband and intraband 
contributions 

, . inter , intra /q (:\ 

Xu± - Xu± + Xu± (A°J 

where we found that X l "± er ~ 0(N 2 ) and Xu™l™ ~ const + O(Nq). So \u± is finite at the transition, as it should be 
when reaching a Fermi liquid. However, the assertion by Millis j|] that the spin-wave velocity should be given by a 
hydrodynamic expression p7J c 2 w — p s /xu± is seen to be incorrect: c sw = c is in fact finite at the transition. The 
correct relationship actually turns out to be c 2 w — p s /x l u± r ■ 



B. Type B 

The non-commutativity of the limits q, lo — > and Nq — ► now leads to much more complex behavior than for 
type A. Consider first, the transverse, staggered susceptibility at q = 0; in this case a relatively complete evaluation 
is possible, and we find 



/ ^ No 

x«±(o,«) ^ 

PS 



vW 



+ 0(w 2 ) 



(3.7) 



where p s ~ Nq ~ (<? c — <j>). The first, singular term arises from interband processes and is due to the integration 
over the region of momentum space where the Fermi surface crosses the magnetic Brillouin zone boundary; there is 
no such point for type A systems. The second term is the regular contribution from the reminder of the momentum 
space, it is similar to that found in type A systems. It is clear that the behavior of \s± is dramatically different for 
u> < No and u > N . 

The q dependence of Xs± only has a relatively simple form in the limiting regions. For q, to <C Nq we have 



Xs±(q,u) 



Nl 

Ps 



N Q 



{biu - ib2ip(q,uj)) 



(3.8) 



where c, b± and bi are constants, and ip{q,uS) = ip(\q x + Qy\, \Qx ~ 9y|i w ) is a symmetrical function of its first two 
arguments, which at q ^> ui/t is linear in q. The explicit form of ip is presented in the appendix. We see that the 
the lowest-energy excitations are overdamped even in the magnetically ordered phase, although we indeed have a 
Goldstone mode at q = 0. Also, the renormalized spin-wave velocity (defined in the limit of small 62) decreases as 



N 



1/2 



In the opposite limit, q,u> ^> Nq (and for all q,u at the transition point) both terms in the last bracket 



in (3.8) become imaginary, and we find 



Xsx(q,v) ^ (q 2 ~u 2 /c 2 -iwy)' 



(3.9) 



The magnon excitations are now clearly overdamped. In the SDW approximation we found that 7 does not depend 
on the ratio w/q. In the quantum disordered phase we then have 



Xs 



(q 2 -L0 2 /c 2 + iu; 7 + A 2 /c 2 y 



(3.10) 



where A is now a pseudo-gap and vanishes at the transition point. Notice also that c is not the velocity of the longest 
wavelength spin-wave excitations in the ordered phase, though indeed it has the same order. 

The same form of susceptibility was proposed earlier by Barzykin et al. M on phenomenological grounds. We will 
use Eqn. (3.10) below to motivate the field-theory for the critical point and the paramagnetic phase. 

As for type A, the transverse, uniform susceptibility is found to have two contributions from interband and intraband 
processes which now take the form 



,, inter , intra 

Xu± - Xu± + Xu± 



(3.11) 



where now X^ er 



O(Nq) and X % uT a ~ const + 0(Nq). Again, x«± is finite at the transition, and the relationship || 



Ps/Xu± does not work as it would predict c s 



Nq. 



practice, however, it is difficult to find conditions when one can neglect the 62 term in (i 



The correct relationship is cj w -- p s /x"i er - ^ n 

Magnon excitations 



are then overdamped, and the issue of the behavior of the spin-wave velocity near the transition is not that relevant. 
Notice also that as shown in the appendix, the O(Nq) contributions in x^"f er and xL"l™ cance l each other, so that 
Xu±= const + 0(N 2 ). 



C. Relationship to results in the Shraiman-Siggia model 



The main weakness of the above mean-field SDW computations is that they underestimate the contribution of 
magnetic fluctuations, which may destroy long-range order at a smaller doping. Such fluctuations may be more 
completely accounted for by studying the phase transition in the Shraiman-Siggia (SS) model 0, which however has 
other weaknesses to be described below. In this subsection we will review the results of such a study JT3] , compare 



them with the above SDW results, and make some speculations on how they may be reconciled. Readers not interested 
in this issue can skip this subsection without any loss of continuity. 

The SS model is a continuum theory of interactions between mobile holes and spin fluctuations in a t-J model 
i.e. a model in which strong local repulsion between electrons allows one to project out states with more than a 
single-electron on a site. In the language of the Hubbard model, there is a significant band-gap between the lower 
and upper Hubbard bands, and the SS model focuses exclusively on physics within the lower band. The lower and 
upper Hubbard bands can be identified with the valence and conduction bands, respectively, of the SDW theory 
discussed above. Furthermore, in the SDW mean-field theory, the minimum, direct band-gap between the conduction 
and valence bands is proportional to the Neel order parameter No, and therefore vanishes when long-range order 
disappears. This is, however, not the case in the SS model, in which No and the band-gap are independent parameters, 
and the band-gap is assumed to remain finite at the point where No vanishes. This is a key difference between the 
two approaches. 

Let us now approach the magnetic transition from the ordered side in the SS model. A simple computation, similar 



to those in Ref 12 1, shows that the transverse, staggered susceptibility has precisely the form ( |3.5| ) obtained above 
in the SDW theory of a type A transition; in particular the r esu lt Im(x s j_) _1 ~ NqUQ is also obtained in the SS 



model (thus the proper interpretation of the factor of Nq in (3.5) in the SDW theory is that of the square of the 
magnetization order parameter, rather than the valence/conduction band gap). At the point where iVo vanishes, 
computations in the SS model show that Q Xs± ~ {q 2 — lu 2 /c 2 + ia'uj 2 q)^ 1 . In the mean-field SDW theory above we 
have Xs_l ~ {q 2 — uj 2 /c 2 )~ l with no damping term; however it is quite reasonable to expect that including higher-order 
paramagnon fluctuation corrections in the SDW theory will lead to a damping term rather similar to that obtained 
in the SS model. So far, therefore, the spin-fluctuation properties of the magnetic transition in the SS model are 
essentially identical to those obtained in the SDW theory of the z = 1 transition in type A systems. 

However, differences do appear when we consider the fermionic excitations. The type A transition in the SDW 
theory has a large electron Fermi surface wholly within the magnetic Brillouin zone, which changes little between the 
two phases on either side of the critical point. On the contrary, the Fermi surface of the SS model, on the ordered side 
of the critical point, consists of elliptical hole pockets at the boundary of the magnetic Brillouin zone. The fate of the 
elliptical hole pockets in the quantum disordered phase of the SS model is not at all clear |13] . Below we present a 
reasonable, but speculative, scenario: we propose that in this phase, the Fermi surface is large (i.e. encloses a volume 
equal to the total number of electrons), as it was in both types of SDW transitions. However the quasiparticle residue 
is very anisotropic so that at the transition to Neel state, the residue vanishes everywhere except for the regions which 
surround the hole pockets in the ordered phase [E8[ . This implies that the critical theory of the z = 1 transition in the 
SS model in Ref |12| remains essentially correct. There may also be states in the Hubbard gap with their residue also 
vanishing at the transition to Neel order. If this large Fermi surface intersects the magnetic Brillouin zone boundary 
in the disordered phase, then there will be a finite T = damping, 7, at q = in the quantum-disordered phase; 
however 7 will vanish faster than the gap, A, as one approaches the critical point, thus behaving as a dangerously 
irrelevant perturbation (i.e., perturbation which is irrelevant very near transition but leads to a new physics at some 
distance away from the transition). Alternatively, the large Fermi surface with anisotropic quasiparticle residue may 
be entirely within the magnetic Brillouin zone; in this case the properties of the quantum disordered phase will be 
very similar to those of a type A transition in SDW theory. 

IV. MODEL FIELD THEORY FOR TYPE B 

This section will present explicit computations of the scaling functions of Sec O in a model field theory appropriate 
for a mean-field type B transition. The motivation for the model follows the logic of Hertz ]15|: the action, S, for 
paramagnon fluctuations in the quantum disordered phase should have a propagator which reproduces the mean-field, 



type B, spin susceptibility in (B.E). The partition function of S in Matsubara imaginary time then has a form 



Z= Vn{x,T) 5(n\x,T)-l) exp(-5[n(a;,r)]) 
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S:= ^ T E/^|l^,^)| 2 ( 9 2 +^/ C g+7KI)- (4.1) 



The action is written in Fourier space, and ui n is a Matsubara frequency. The vector field n(x, r) represents the local 
orientation of the antiferromagnetic order parameter; we will allow n to have N components to allow a subsequent 
large N calculation. We have chosen to implement a fixed-length constraint ft 2 = 1 to mimic interactions between 



the paramagnon modes. This restriction is however not crucial and identical scaling results would be obtained in 
a model with a more conventional (u/2N)(n 2 ) 2 interaction. The scaling results appear a little more directly in the 
fixed-length model. The first two terms in S are those found in the usual O(N) non-linear sigma model. The last 
term arises from the damping induced by the fermion particle- hole pairs. The action S was first explicitly written 
down by Liu and Su fa] — they used a two-component microscopic model in which the spins and fermions are locally 
independent degrees of freedom, and then integrated the fermions out. 

We expect S to provide a reasonable description of the temperature dependent crossovers in the quantum disordered 
phase. However, close enough to the critical point in the quantum disordered phase, the validity of S appears to us 
to not have been established convincingly. It is likely that it will be necessary to account for fluctuations involving 
incipient formation of the spin-density-wave gap over portions of the Fermi surface — note that these gaps form at 
precisely the same points on the Fermi surface which are responsible for the low frequency damping on the disordered 
side; these issues will be examined in more detail in subsequent work. It is also easy to see explicitly that S certainly 
breaks down on the magnetically ordered side of the transition. On this side, the spin fluctu atio n propagator changes 



substantially when q and w become smaller than the order parameter Nq (see e.g., Eq. (3/7)), and these changes 
cannot be implemented by simply introducing a condensate of the order parameter into S. 



In the remainder of Sec. [V we will focus exclusively on the model S. We will describe its asymptotic critical 
behavior, and compute its scaling functions in a large N expansion. 

It is also interesting to note in passing that the \u>\ dissipation in S is similar to that used in the macroscopic 
quantum tunneling literature |29(| . There the emphasis is on the consequence of this dissipation on the quantum 
mechanics of a single, heavy particle, whereas we are studying its effects on a field theory describing a large number 
of degrees of freedom. 

Let us first discuss some general scaling properties of S from the vantage point of the primary, z = \ fixed point. 
This fixed point clearly lies on the submanifold of the parameter space with 7 = 0. The term proportional to 7 is 
clearly a relevant perturbation at the primary fixed point. Moreover, \w n \ is non-analytic in frequency and leads to 
a long-range — 1/t 2 interaction in spacetime. Theories with long-range interactions of this type are familiar in the 
context of finite temperature classical phase transitions [ p0[ and many results can be transferred over. In particular, 
the \ui n \ term gets no singular loop renormalizations, and the scaling dimension of 7 is exactly 1 — 77. Therefore, at 
T = 0, the response functions will therefore be scaling functions of the ratio 7/(5 — <? c )^ 1_r " I/ ■ The fully renormalized 
parameters introduced earlier will have a different dependence on the bare coupling constants depending on the value 
of this ratio; this dependence can be deduced by standard scaling arguments. For 7 <C (g — g^ 1 ~ T, ' v we will have 
dependencies characteristic of the primary fixed point for which 

Z~ constant, A~{g-g o y, T ~ 7(3 - g e )"" '. (4.2) 

In the opposite limit 7 ^> (g — g c ) we are much closer to the critical point and dominated by the secondary 

behavior; in this case we have 

Z~ constant, A ~ (g - jji/ywW/P-^ r „ y/d-j). (4.3) 

Explicit results for scaling functions of S can be obtained in a 1/iV expansion. The methods are very similar to 
those discussed at length in Ref. Q for 7 = 0; one only has to add a 7|w n | to each propagator. We will therefore 
merely present the final results. At N — 00 we found 

$f =oc (q, 57, A, T) = -5 5 J =- . (4.4) 

s W ' ' ' ' q 2 -ZJ 2 -tTZJ + m 2 (A 1 T) K ' 

We chose the definition of the parameters earlier so that as T — ► 0, m = A. For finite T, m is given implicitly as a 
universal function of A and Y by 




_ f m \ - /A \ [°° dtt , / TO 



n 2 



where the value of the arctangent runs from to 7r as f2 runs from to 00, and 

{Ax 2 - l) 1 / 2 tan- 1 (4a; 2 - l) 1 / 2 for x > 1/2 

<Kx) = { (1-4^1/2 1 _ fl _ 4 2U/2 • (4-6) 

IL -T J -'<* U?i-tr l * f0 " sl/2 



Despite appearances to the contrary, <f>(x) is analytic for all real x, including x — 1/2. In the limit r — » 0, ( |4.5| ) 
reproduces the z = 1 result of Ref M: m = 2sinh~ 1 (e A / 2 /2). In the opposite limit r — > oo we expect properties of 






the se con dary fixed point — we see from (2.5) that m 2 /T should be a function of A /T. This is exactly what we find 



from (4.5) which reduces to 



m 2 -A 2 m 2 T 2 A 2 , T 2 n f°° tan' 1 (SI/ (m 2 /T)) 

= + ^log— -^lo g= 2 =2 / dSl \ n ' >> 4.7 

r r m l r a h e - i 

for large T. The only violation of the A /T scaling in the above comes from the logarithms, which was also expected; 



the arguments of the logarithms however remain universal. Analysis of the solution of (4.7) for m can be shown 
to yield the same crossovers (including and upto log log corrections) as those discussed in Ref. la] in the quantum 
disordered region for d = 2, z = 2. 

The calculations of the 1/N corrections also follow Ref j|. The evaluation requires substantial numerical com- 
putations which are currently being carried out. However, we also know from Ref 0j that at least at T = 0, these 
corrections were quite small, and that the N — oo results were satisfactory for most purposes (a notable exception is 
X'''{u) for r -C A, for which 1/N corrections are relevant). The agreement of the large r, N — oo results with those 
of Ref. H is also encouraging. 

A number of experimentally measurable quantities can be deduced from the above results. For example, the 
antiferromagnetic correlation length is 

and the NMR relaxation rates (for a review, see, e.g., pW), upto known prefactors associated with hyperfine coupling 
constants, are |7|,|32|] 



i r 


1 1 


T ^ 2 


T 2G ^ mT 



(4.9) 

We plot in Fig [| the universal crossover function for hc/t^A as a function of fc^T/A for various values of T/A. For 
f « A we have two main regimes of temperature: 



(a)/c B T>A C = 21og 




«,)k B T«A i-?|l-l=fH-A)t;^J)). (4.10) 

Regime (a) is the z — 1 quantum-critical behavior (lj,y], and (b) is the z — 1 quantum-disordered. At very low 
temperatures, the damping gives rise to a power-law rather than exponential behavior of the correlation length; these 
power-laws are the same as for quantum disordered, z = 2 behavior, discussed below. However, unlike for "pure" 
z = 2 relaxational behavior, x"(w) has a knee at u> ~ A ^> T. Based on the presence of such knee (shown in Fig.|l|b), 
we identify this regime as quantum disordered, z = 1. 
For r>A there are three sub-regimes: 

(a') fc s T»r e = 21og' V r>+ ' 

(b') ^r « k B T « T e = /: 

(c>)k B T«¥- £ = ^ [i-____ (1^) (A.)]. (4.11) 

with fi(x) a very slowly (logarithmically) varying, numerically calculable function of order unity. Now regime 
(a') is z = 1 quantum-critical, (&') is z — 2 quantum-critical, and (c') is z = 2 quantum-disordered. In 
the negligibly small subregime of (&') where \og\og(T /k B T) ^> 1, while maintaining k B T ^> A 2 /r we have 
£ = ?7c/(7rrfc s T) 1 / 2 [log(r/fcBT)/loglog(r/fcBT)] 1 / 2 . Most of the above asymptotic results are not terribly useful 
at any realistic value of T/A, and the exact values plotted in Fig H should be used for experimental comparisons. 
The above regimes are sketched in Fig Bl which will be discussed in Sec [v| 
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A. Uniform susceptibility 

Consider now the response of the nearly antiferromagnetic Fermi liquid to a uniform magnetic field. Unlike the 
staggered susceptibility which diverges at the transition, the uniform susceptibility, \ Ul is expected to remain finite. 
As a result, the crossover behavior of \u is more complicated and corrections to the leading scaling result now play a 
significant role. 

The following results are for type B systems; however they can also be applied to type A systems after imposing 
7 = T = 0. 

We need the modifications to the action S, in the presence of an external field H . These corrections are equivalent to 
evaluating the effective 3-point and 4-point couplings between the antiferromagnetic order parameter and the uniform 
magnetization after integrating out the fermions; however, as we are mainly interested in the case of a uniform field, 



we will just carry the exact field dependence through in all the fermion propagators (see Sec A3). 

In insulating antiferromagnets, the coupling to the field was given exactly by a symmetry-related argument |||] 
with no additional coupling constants appearing. These symmetry arguments must now be applied with more caution 
because of the presence of gapless fermionic excitations in the quantum-disordered phase. In particular, we need to 
account for the continuous change in the Fermi liquid ground state in the presence of H associated with the expanding 
(shrinking) of the up-(down-)spin Fermi surface. 

Consider, first, the effect of H on the \u> n \ term in S. At small oj n , this term arises from the absorption due to 
particle-hole pairs with momentum Q, made up of electrons and holes just above and below the Fermi surface. A 
field, H, will shift the bottom of the band of the up and down spins in opposite directions. In a model with a flat 
density of states, such a shift should have no effect on the particle-hole spectrum. So for small H , there is no change 
in the \u) n \ term. At larger H there will be a change, because there is always some structure in the density of states: 
this will show up as a field dependence in the value of the coupling 7: 7 — ► 7(1 + XH 2 ) where A is a small coupling. 
This is verified in an explicit computation of the field dependence of S discussed in the appendix. 

The computation in the appendix also finds a precession term: 



1. ,. v / j2 . » - dn 



—zCl + ajdrxdTH-rix—, (4.12) 

get J or 

where we have absorbed in the definition of H a factor gcfJ-B/^ (gc^B is the gyromagnetic factor). The prefactor of 



(4.12) has been written such that the coupling a = in the insulating limit; in this limit the precession rate of the 



spins is known exactly and hence (4.12) has no new coupling constants. In the doped system, there is a correction to 
the precession rate from internal fields generated by the polarization of the fermions, and this is represented by the 
new coupling a. If the system is not too strongly doped, we can expect that a is not too large. 

Finally, there is a term (also derived in the appendix), ~ (ft x H) 2 , which imposes an energetic preference for the 
relative orientation of the antiferromagnetic order parameter and the magnetic field. 

All three terms discussed above induce non-universal corrections in the T = value of Xu- These may be interpreted 
as corrections to the fermionic Pauli susceptibility from higher-order interactions between the paramagnon modes. 



Things do simplify, however, when we consider the T dependent part of Xu- Then the precession coupling in (4.1 2j) 



turns out to be the most important in many cases. The T-dependent corrections from the H 2 dependence of 7 can be 
shown by standard scaling arguments to be subdominant near the pri mary fixed point; near the secondary fixed point 
this term turns out to have the same T dependence as that due to ( f.I2| ), but with a prefactor of (7/i) 2 . Finally, 
the last (n x H) 2 term yields no temperature dependence in a model with a fixed- length constraint. In a soft-spin 
theory, scaling arguments show that this term is also subdominant near the primary fixed point; near the secondary 
fixed point, and going beyond N = 00 limit for soft-length spins, one obtains a T dependence similar to that due to 
the H 2 dependence of 7. 

We see therefore, that if 7 is sufficiently large, all three terms contribute roughly equally to the uniform susceptibility. 
The temperature dependence of Xu is then non-universal, and reliable predictions for experiments are difficult. To 
keep our discussion simple, and also because we think that this is most realistic experimentally, we will only present 
explicit results for x« for the case when 7 is not too large. We will therefore only consider the T dependence of Xu 



induced by the term ( 4.12 ) which domi nates near the primary fixed point 



We now consider the consequence of (4.12) on the T dependence of Xu- Application of our earlier scaling arguments 



and the results of Ref |4j allow us to deduce the following crossover scaling function for Xu ■ 

Xu(T) - Xu(T = 0) = (f + c/)M^„(A,r) (4.13) 



ff 



where a' is a non-universal constant, which vanishes when a = 0. The function $ u is a completely universal scaling 
function which we will compute below at N = oo. The only non-universalities of the above result are therefore related 
to the coupling a', and the value of Xu(T = 0) which is dominated by the Pauli susceptibility of the fermions. 
The N — oo computation of <I> U for S can be performed as in Ref M; we find: 



<=°°(A,r) = ^ 




2A 2 -r 2 . /A 



4A 



fir 



ft 3 



oil 



-1 



T^2, 



(m 2 -ft 2 ) 2 + r ft 2 



(4.14) 



where m has to determined from (4.5) as a function of A and T, and the function <p(x) was defined in (4.6). A plot of 
the universal contribution from $ u to \u is shown in Fig |5J, with the non-universal prefactor 1 + a' dropped. There 
are several different regimes, similar to those found in the correlation length. Dropping the 1 + a', we state the leading 
behavior of Xu(T) in these regimes. For T <C A we have the two regimes of temperature: 

(a) k B T » A XuiT) - Xa (0) = ^ log ( ^±l) **£ 



(6)fc B T«A Xu(T) - xu(0) = A (cxp 



A \ 7T r /fc R r 



/crT/ 6 a v a 



(4.15) 



For r 3> A the three sub-regimes are: 



(a') k B T > T 



m m y/E VE+l \ k B T 

XuiT) - x„(0) = — log — - — — r - 



A 2 

{b') — « fc s r « r 



(c') fc B T « 



A 2 



Xu(T) - Xu(0) = f 2 



Xu(T)- X u(0) 



A:rT\ /crT 



r J c 2 



fc B T 



61og(r/A) V Ac 



(4.16) 



with fi a function similar to the function /i above in (4.11). The regimes (a), (6), (a'), (6'), (c') are the same as 
those described for the correlation length (notice that in the regime (b), the damping term gives rise to a power-law 
behavior at very low T). Again, as with the correlation length, most of the the asymptotic expressions are not useful 
for realistic T/A and numerically determinable values of the crossover function should be used. Within the regime 
(&'), it is possible, as was the case with £, to have a negligibly small sub-regime where loglog(r/fcsT) ^> 1 while 
k B T » A 2 /r; here we find X u{T) - % u (0) = (fc B T/ 7 rc 2 )(loglog(r/fc i3 r))/log(r/fc i 3r). 



B. Specific Heat 

The leading contribution to the free energy density of S is strongly divergent in the ultraviolet ~ A 3 , where A is 
an upper cutoff in momentum space. In the model with 7 = 0, it was found in Ref. H that a single subtraction of 
the T = value of the free energy density was sufficient to make the T dependent remainder finite. There were no 
terms diverging with powers of A smaller than 3 — a direct consequence of the hyperscaling properties of the original 
theory. For finite 7, we find here that this simple property does not hold: a single subtraction still leaves a term that 
is only logarithmically dependent on the upper cutoff. This is perhaps not surprising as the secondary fixed point 
violates hyperscaling. While most of this violation is actually cutoff by the crossover to the primary fixed point at 
high energies, some manages to survive. 

We will restrict ourselves here to presenting the T dependence of the free energy density, J 7 , of the action S in a 
N = 00 calculation. We found 



^(T)=^(o) + ^p%>(A,r) 

where the dimensionless function 'J, instead of being fully universal, has a logarithmic dependence on A/T: 



(4.17) 
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(4.18) 



where to is determined from (4.5) as a function of A, T and the function <fi is defined in (4.6). 

The asymptotic limits of the contribution of \& to the specific heat Cy = —Td 2 TjdT 2 were determined in a manner 
similar to £ and Xu- For T <IC A we have: 



(a) fc B T > A 

(6) fc s T < A 

For r ^S> A the three sub-regimes are: 

(a') fc s T > r 
A 2 

(b') — « k B T « r 

A 2 
(c') fc s T « — 
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(4.19) 
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(4.20) 



with /3 a function similar to the function /i above in (4.11). The regimes (a), (6), (a'), (6'), (c') are the 
same as those described for the correlation length. Within the regime (&') above, it is possible, as was the 
case with £, to have a negligibly small sub-regime where loglog(r/fcsT) 3> 1 while k B T ^> A 2 /I"; here we find 
C v = (N/12){rk B T/(hc) 2 )\og((hcA) 2 /k B TT). 



V. CONCLUSIONS 

We begin this concluding section with a simple, qualitative discussion of the physical picture behind the compu- 
tations in this and related, previous, works on the spin fluctuation properties of the not too strongly doped cuprate 
compounds. A discussion of the relationship of our ideas to experimental systems will follow. Finally we will discuss 
some open theoretical issues. 

It is very useful to think in terms of the physics at different length and energy scales. On the whole, we may assume 
that the relevant energy scales decrease uniformly with increasing length scales, so the two can be used interchangeably 
to move between the different regimes. In what follows, we describe the sequence of crossovers as one moves from 
larger to smaller energy scales, or equivalently from smaller to larger length scales. 

At the very largest energy scales, the behavior is dominated by lattice scale physics, which is inherently non- 
universal. At slightly smaller energies, provided the doping is not too large, one can neglect the effect of mobile 
carriers. It was argued in Ref that spin fluctuations at these scales are "quantum- critical" and well described 
by properties of the z — 1 critical point in the 0(3) non-linear sigma model which separates the Neel ordered and 
magnetically disordered phases. This proposal has been the source of some controversy in the literature, although 
some fairly convincing evidence has appeared in recent high temperature series studies [p4|. The excitations in this 
2 = 1 quantum-critical regime are neither spin waves of the ordered state, nor 5 = 1 magnons of the disordered state, 
but form a critical continuum. 

The subsequent crossovers at smaller energies depend on the doping concentration. At very low doping, the presence 
of T — long-range order in the Heisenberg antiferromagnet becomes apparent; there is then a crossover from z = 1 
critical fluctuations at larger energies to the Goldstone spin-wave modes of the ordered state at lower energies. At a 
slightly larger doping, the T = long-range order is destroyed. The crossover at smaller energy scales is then into 
the quantum-disordered state of the sigma model, where the excitations are then gapped, triply degenerate, spin-1 
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magnons. At this same doping, and provided the geometry of the Fermi surface is of type B, there is a last crossover, 
at an even smaller energies, to a state where the mobile carriers induce an important T = damping of the magnon 
excitations p5|,p[; this damping leads to to sub-gap absorption and power-law (in T) behavior of observables. All 
power-laws are equivalent to those of quantum disordered, z = 2 regime, discussed below. 

This sequence of crossovers has been described for the case in which the damping is not too large: T < A in the 
notation of Sec O. 

At even larger dopings we may have T > A; then the intermediate quantum disordered, z = 1, regime of a pure 
sigma-model is not realized, and the system undergoes crossovers from z = 1 critical regime to z — 2 critical regime, 
and then at even smaller energies to z = 2 quantum disordered regime with relaxational behavior of spin excitations. 

The scaling regimes and associated crossovers are sketched in Fig pi We emphasize that all phase boundaries 
between phases on this sketch are smooth crossovers, and their precise positions are therefore not meaningful, while 
the sequence of crossovers and overall structure of the diagram is meaningful. In terms of this phase diagram, the 
Shraiman-Siggia model (type A), corresponds to r/A = |4|Jl2|] (modulo dangerously irrelevant damping terms); 
models with primarily relaxational behavior of spin excitations (type -B), similar to that observed in Y \Ba2Cu3O7 , 
to r/A > 1 ifulpl; and, finally, models with moderate T/ Delta < 1 exhibit spin pseudogap behavior similar to 
that observed in the underdoped Y-based materials, such as Y Ba2Cu30e.63 and Y BaiCunQz fg,0]. The possible 
relationship of the phase boundaries to doping levels in the cuprate materials will be discussed in more detail in a 
subsequent publication. 

Going to a finite T introduces an additional degree of complication which substantially changes the crossovers 
described above. However, all of these T dependent crossovers are contained in the universal crossover function <& s in 



Eqn (2J3), which was computed in a large N expansion in this paper. The main results can be understood by keeping 
a simple rule of thumb in mind: the primary effect of a finite T is to thermally quench the excitations at the energy 
scale ksT, so that the crossovers at scales below fcgT no longer occur. It is therefore possible to make experimental 
predictions for temperature dependences that are remarkably universal [[|-|| . 

We now briefly describe the relationship of our results to experiments on high-T c oxides. We leave detailed quan- 
titative comparisons for a separate publication which will also contain computations of \/N corrections; here we 
underline only the main qualitative points. 

La2- x Sr x CuO±: The Fermi-surface of La-based materials has not been measured experimentally. Apparently, they 



are close to type A 1 20 in the terminology of Sec El and we expect to see only the features of the z = 1 behavior. 
However, these materials also show nontrivial spin correlations in the metallic phase which will certainly effect the spin 
fluctuations at low enough temperatures |3q] . The primary effect of mobile carriers on the uniform susceptibility should 



be a temperature- independent additive contribution, which increases with doping. We have suggested earlier plp2 36 
that this behavior may have been observed at a few percent Sr concentration. 

YBa2Cu^O-j: Photoemission experiments |2J| indicate that the Fermi surface is large and belongs to type B. 
Both our and earlier |16|] analysis of the NMR data is roughly consistent with type B scaling behavior at A/T ~ 1: 
the measured \jT\T monotonically decreases with T, and the ratio T\TjT^ is temperature independent p7|| . A 
quantitative comparison with the m e asure d uniform susceptibility |3q| is difficult because nonuniversal contributions 
to susceptibility, neglected in ( 4.15 ), J4.16| ), are relevant for this fully doped material 



YBCL2CU3O6.6 and YBa2Cu4 : Os: For both of these materials, the measured Fermi surface is large and belongs to 
type B ]3^,^(J. The measured \jT\ |4l],|^] and uniform susceptibility [g8| decrease rapidly as temperature decreases 
below T ~ 150.K\ Such behavior is characteristic of the quantum-disordered z = 1 regime and we therefore expect 
that for these materials, r/A <C 1. Note however, that we have shown explicitly that due to the finite damping, this 
decrease, particularly for the uniform susceptibility, is not purely exponential as it is in phenomenological "spin-gap" 
models. We also note that it is important to include 1/N corrections Q in describing the behavior of \jT\T in this 
regime. At high temperatures, this system is expected to crossover to z = 1 quantum-critical behavior, where T\T 
is linear in T and T1T/T2 — const HJq]; this behavior appears to have been observed in both YBa2Cu-iO§.§ and 
YBa 2 Cu 4 8 H||. 

We conclude this paper by recalling and raising some open theoretical questions. The explicit computation of 
temperature-dependent crossovers has been restricted to the quantum disordered phase, and not too close to the 
quantum-critical point. We suggested that S may fail close enough to the critical point, and saw explicitly that it 
failed badly in describing the low energy excitations on the ordered side. A drastic change in the action seems to be 
necessary; probably, much more drastic than adding some additional dangerously irrelevant coupling which becomes 
important only the ordered side. In fact, any procedure which uses a truncated power series expansion in ft for the 
action appears doomed: no such expansion will correctly capture the appearance of gaps on portions of the Fermi 
surface. It seems that the only natural way of describing the crossover on the ordered side is to work in a theory in 
which the fermions are not completely integrated out. It is then possible to obtain an action which has only regular 
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terms, and which holds on both sides of the transition. A possible form of the action involves two species of fermions 
(corresponding to pairs of points on the Fermi surface separated by Q) coupled to the n-held. Computing crossovers 
for such an action remains an important open problem. 

A complementary open problem applies to the analysis using the Shraiman-Siggia J12[ model. In this case the 
ordered state and its crossovers are reasonably well studied. However it is also clear that the model cannot apply in 
detail on the quantum disordered side, and the crossovers between the ordered and disordered sides are not completely 
understood. 

Finally, we make a few remarks about the effects of disorder. Hertz |15 did begin to address the consequences of 
disorder, but his analysis is seriously incomplete in all cases. One sign of this is that his value of v (which equals 
1/2 in all the cases considered by him) violates the inequality v > 2/d [H where d < 4 is the spatial dimension. 
It is easy to determine the source of this violation: the most important perturbation on a soft-spin version of an 
action like S is a random- "mass" term J d d xdr ~ m 2 (x)ft 2 (x , t) which accounts for fluctuations in the local value 
of the critical coupling at which the magnetic order is destroyed. It can be easily shown that this term is relevant 
about the pure fixed point as long as v < 2/d. A study of S, with an additional random mass, along the lines of the 
work of Boyanovsky and Cardy |4f| should be quite straightforward. However this approach involves the potentially 
dangerous procedure of expanding about a problem with e T time dimensions, and a method which works directly with 
e T = 1 would be preferable. 

Note added: Just before this paper was submitted, we received a preprint from L. B. Ioffe and A. J. Millis |4q | 
which addressed the T dependence of the uniform susceptibility near the type B fixed point. Their expression for Xu 
in terms of correlators of n is in general agreement with our discussion in Sec IVA. However, they focused on the 
temperature dependence of Xu at very low T (logT/T ^> 1) near the type B fixed point; we chose to focus on the 
term important near the primary type A fixed point — the latter term gives a universal (up to an overall prefactor) 
contribution to Xu which should be dominant at all T unless T/A is very large. 
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APPENDIX: SUSCEPTIBILITIES IN SDW THEORY AT T=0 

1. Staggered Susceptibility 

In the SDW theory, the spin susceptibility is given by a ladder series of bubble diagrams. At T — 0, one fermion 
in the bubble should be above the Fermi surface, and one below. In principle, the solution of the ladder series in 
the magnetically ordered state is a 2 x 2 matrix problem, as one has to consider bubbles with momentum transfer 
and Q Hffl. For our considerations, it turns out that all terms with momentum transfer Q disappear at the 
critical point, and we have checked that they only account for small corrections to the expressions obtained below. 
We therefore neglect terms Xo (<iS q + Q,u>) with the momentum transfer Q in which case the total transverse dynamic 
susceptibility is given by 



x + (g,w) 
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Here and below we set h = 1 . We will now study the behavior of this result near the magnetic transition where Nq 
becomes small. 
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Near q = Q, the coherence factors can be simplified to 

N 2 

N 2 
(H+e- k+q f^y^O((Q- q f). (A3) 

We first consider the behavior of dynamic susceptibility at the antiferromagnetic momentum, q — Q. Here e k + e fe+„ = 
0, E k + E k+ = 2E k , and one obtains using the self-consistency condition 

xt(Q,«») = rr + JfY,\< - n t) ( ^ + ^— - ^=) ■ (A4) 

[/ N^ \2E k ~uj 2E k +u E k J 
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Expanding this expression in ui and substituting into ( |Al[ ) we obtain 

(y + (Q,u>)) = > — -. (A5) 

The difference between type A and B transitions now becomes transparent. In the first case, the integration over k 
is restricted to a region which is located entirely inside the magnetic Brillouin zone. Accordingly, E k remains finite 
when No tends to zero, and at the critical point we have (x H (Q,o;)) _1 ~ oj 2 . For type B, however, the allowed 
region of momentum integration includes the vicinity of (xq,^ — Xq) and symmetry related points where the Fermi 



surface crosses the magnetic Brillouin zone boundary. At each of these points, E k — No, i.e, the denominator in (A.5) 
diverges as Nq — > 0. Expanding E7 near these points and performing the momentum integration, we obtain after 
some algebra 

(X+ - (M) - I= _„£(_!= 4) (A6) 

where D = t/87r|i'sinxo sin2xo|, and the uj 2 term comes from the integration over the regions far from (xo,7r — xq). 
We see that the leading term in the expansion in w is now u> 2 /No- As one approaches the critical point, the first term 



in (A6) becomes purely imaginary, and at Nq = one has 

(xMQ^r 1 = -D^r (iu> + ^) • (A7) 

We now consider the form of the susceptibility at finite q which we will understand as a deviation from Q. The 
calculations here are straightforward but tedious, so we skip some of the details. First, we found that at w = 0, 
all potential terms of the form q 2 /No cancel each other for both type A and B— type transitions. The expansion of 
(X"' (9;0)) _1 over q is therefore regular for both cases. Second, we found that when both to and q are finite, the 



second piece in (A2) has an imaginary part which for type A transition behaves as 



Im( X +-)- 1 ^^iV 2 . (A8) 

We see that the damping of magnetic excitations disappears as No — > 0. 

For type B transitions, imaginary part of the susceptibility is a function of q/No and to /No- For q, lo <C No, we 
obtained 

U 2 
Im(x + ) 1 =D — il>{\q x +q y \,\q x -q y \,u) (A9) 

where D = £/[87r|£'sin2xo|], and ip looks particularly simple at t' <C t 

where Co = u>/(2t sinxo). We also obtained expression for ip for arbitrary if /t, but it is more cumbersome and we refrain 
from presenting it here. In the opposite limit, q,u> ~^> No, or equivalently when No vanishes, the q/No dependence 

1G 



in (A9) transforms into a constant term, and both pieces in (A2) contribute linear in u> terms into the imaginary 
part of x H • The numerical factors for each of the two pieces contain combinations of 9— functions of the form, e.g., 
9(At sin xo\q x + q y \ — \u> + 4|i' s'm2xo\(q x — q y )\). However, we have checked explicitly that the sum of the two terms 
does not contain 6 functions. Specifically, we obtained at Nq = and arbitrary ratio ui/q 



Im( X +-)- 1 = J D 



El 
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Combining the above results, one obtains the expressions presented in Eqns (pjj), (|3.8|), (p.9| 



2. Uniform susceptibility 

The uniform transverse spin susceptibility is defined as 

\2 u„, ,.+ 
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(A12) 



where qg^b is the gyromagnetic factor and a 2 is the volume per spin. In the q — ► limit, the coherence factors are 
simplified and we obtain x± = (.9GA i s/a) 2 Xo/[2(l — f Xo)] where 



*°-n2^ 



E7 



El 



^o AT 4^ R- 



lim 



'fc+9 



jv'd 

h k+q 



l k+q 



E 



k+q 



El 



(A13) 



As one approaches the critical point, the first term in ( |A13| ) disappears and x± takes the familiar RPA form of 
magnetic susceptibility in a Fermi-liquid, which is finite at T = 0. For the type A transitions, the first term vanishes 
as Q(N 2 ), a nd for the t ype B transitions it vanishes as 0(Nq). The uniform susceptibility then has forms presented 
in Eqs. J3.6| ) and ( 3.11 ). Notice however, that the full xo (arid hence, x_l) behaves as xo = const + 0(N 2 ) for both 
types of transitions: for the type B transition, linear in N$ contributions from the first and second terms in ( A13 ) 
cancel each other. 



3. Staggered susceptibility in a finite magnetic field 



To obtain the temperature dependent correction to the uniform susceptibility from the fluctuation and interaction 
of the magnon modes, we will need the form of the long-wavelength action for antiferromagnetic fluctuations in the 
presence of the external magnetic field. To derive this form, we need to understand how the staggered susceptibility is 
modified in such a field. In this subsection, we will determine this modification in the SDW formalism for the ordered 
state. We then extend SDW calculations upto a critical point which in turn will allow us to obtain the field-dependent 
part of the action in a disordered phase. 

To be definite, we assume that the field is directed along x axis (antiferromagnetic order parameter is along z axis). 
We also introduce new operators as linear combinations of quasiparticle operators with up and down spins 



Ofc i + au,2 , 

'fe = ^r— ; Wk 



«fc,l — Ofc, 2 



V2 



V2 



(A14) 



(indices 1 and 2 correspond to up and down spin, respectively) . The quadratic part of the Hubbard Hamiltonian now 
takes a form 



E (£fc+Q - H)ct>l +Q 4) k+ Q + (e k + H)iplrp k - N (ct>l +Q ip k + ^{<t>k+Q 



(A15) 



where H = H/2. Clearly, the magnetic field splits conduction and valence bands for up and down spins so we now have 
four different branches of fermionic excitations. However, we see that the pairs of operators (4> k , ipk+q) and (<pk+Q, ipk) 
in ( |A15| ) are decoupled from each other, and therefore the diagonalization is still a 2 x 2 problem. Performing the 
standard manipulations we obtain 
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where we introduced 



Hsvw = Y^ E >W + E iAfa + E Hlk + E 5 k 5lS k (A16) 



E a ' P = 4 ± ^~, i^ = 6+ ± £ fc "+ (A17) 



and E k = a/ A, 2 + (e fe — if) 2 ; E k + = a/A 2 + (e k + H) 2 . The self-consistency condition now takes the form 

! = Ey ' f"fa-"«» + ^-'M . ( A18) 

The computation of susceptibilities proceeds in the same way as without field, the only new feature is the appearance 
of a cross-polarization term Xq z which makes RPA ladder summation a 2 x 2 matrix problem even if we neglect, as 
we did earlier, small terms with a momentum transfer Q. Specifically, we found 
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Below, we restrict the consideration to an analysis of the frequency dependence of the susceptibilities at q — Q. In 
this case, the bare susceptibilities are the following 
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Performing the computations, we found that 

1 oj 2 /.. H_ 

X z v (Q^)^ic 3 u;H/t 3 (A21) 

where c; = Ci(t'/t) and Hq ~ t. For the type A transition, we found that all three coefficients remain finite at 
the critical point. Elementary considerations then show that the dominant field dependence in the denominators 
of (|A2C|) comes from (xo V ) 2 an d is in the form uj 2 H 2 . For the n— field action, this implies that the extra term in 



the presence of H has a form of Eq. (4.12) (to see this, one simply has to compute correlator n z n v ex x zy )- For a 
disordering transition of type -B, the situation is more complex. The most relevant point is that C3 is still finite at 
the transition which implies that there is no term ~ sgn(uj) in the coupling to the field. However, we also found that 



ci and c 2 behave as a = (l/2)D((t/y/4N§ - u> 2 ) + ...), c 2 = D((t/2N Q ) + . . .) where D was defined after flA6| ) and 
dots stand for non-singular (field-dependent) terms. Right at the critical point, c\ — [(i/2)Dt/w +...], and we have 
Xl y ' zz = (1/2U) - iju(l - H 2 /Hl)/(W 2 ) + 0{uj 2 (1 - H 2 /H?)). We see therefore that for B- type transitions, 7 



acquires a field-dependent part: 7 — > 7(1 — H /Hq). Substituting this expression into (A.2C), we find an extra field 
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dependence in the denominator which has the same uj 2 H 2 form as the dependence imposed by (xo V ) 2 , but with extra 
prefactor (7/i) 2 . This implies that for z = 2 transition, there are two different relevant terms in the n— field action 
which describe the coupling of the antiferromagnetic fluctuations to the field. We discuss the relative importance of 



these terms in Sec. IVA. Finally, the effective action indeed contains a term H 2 — (Hn) 2 = (n x H) 2 which can 
be extracted from the longitudinal susceptibility (x xx (0, 0)) _1 ex H 2 . However, because of the constraint, this term 
contributes only to uniform susceptibility at T — which is non-universal and is not considered in this paper. 
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FIG. 1. Sketch of the imaginary part of the local dynamic susceptibility, x'l (obtained by integrating over momenta in the 
vicinity of the antiferromagnetic ordering wavevector) of a nearly antiferromagnetic Fermi liquid in two dimensions at zero 
temperature, (a) The T — damping constant T is exactly 0, as is the case for an insulating antiferromagnet in a quantum 
disordered phase; A is then the spin gap. (b) Shows the consequences of a small value of T; there is now no true gap, only a 
gap- like knee in the spectrum, (c) A large value of V makes the spectrum relatively featureless. 

FIG. 2. (a) Fermi surface of type A: ordering wavevector Q = (n,n) cannot connect two points on the Fermi surface; (b) 
Fermi surface of type B: Q can connect two points on the Fermi surface. Dark regions correspond to occupied electronic states 
in the first Brillouin zone. 

FIG. 3. The evolution of the Fermi surface with doping in the Hubbard model with first and second neighbor hopping. 
Shaded regions correspond to hole pockets and black regions correspond to electron pockets (i.e., doubly occupied states in the 
first Brillouin zone). The description of the figures is given in the text. 

FIG. 4. Universal crossover behavior of the inverse antiferromagnetic correlation length l/£ (measured in units of 1/h c) as 



a function of the temperature T, for a variety of values of T/A. The plot is the TV = oo result contained in Eqns (4.5) and 
([1.8]). There are no arbitrary scale factors, and the scales of both axes are universal and meaningful. 

FIG. 5. As in Fig[|, but now with results for the temperature-depen dent part of the uniform spin susceptibility Xu(T) — X«(0) 



measured in units of (qg Ub/ Hc) . The non-universal constant a' in (1.13) has been set equal to 0. The plot is the N = oo 



result contained in Eqns (4.5) and (4.14). Apart from the non-universal scale factor associated with a' , the scales of both axes 



are universal and meaningful. 

FIG. 6. Scaling regimes are sketched versus the normalized temperature, T/A, and the T — parameter T/A which 
describes type A to type B crossover. While the asymptotic T — > behavior is the same for all nonzero T/A, an important 
distinction between the quantum disordered (QD), z = 1 regime at small F/A and the QD, z — 2 regime at large F/A is the 
presence of a knee (see Fig nib) in x"( w ) at w ~ A in the former case. We emphasize that the T — parameters A and T may 
have a highly nontrivial, possibly non-monotonic, doping dependence, therefore the boundaries between different regimes on 
this phase diagram may have somewhat different structure compared to a phase diagram plotted as a function of doping. 
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